Dilution and resonance enhanced repulsion in non-equilibrium fluctuation forces 
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In equilibrium, forces induced by fluctuations of the electromagnetic field between electrically 
polarizable objects (microscopic or macroscopic) in vacuum are always attractive. The force may, 
however, become repulsive for microscopic particles coupled to thermal baths with different tempera- 
tures. We demonstrate that this non-equilibrium repulsion can be realized also between macroscopic 
objects, as planar slabs, if they are kept at different temperatures. It is shown that repulsion can 
be enhanced by (i) tuning of material resonances in the thermal region, and by (ii) reducing the 
dielectric contrast due to "dilution". This can lead to stable equilibrium positions. We discuss the 
realization of these effects for aerogels, yielding repulsion down to sub-micron distances at realistic 
porosities. 

PACS numbers: 12.20.-m, 44.40.+a, 05.70.Ln 



Forces induced by electromagnetic (EM) field fluctu- 
ations of quantum and thermal origin act virtually be- 
tween all matter that couples to the EM field, since the 
interacting objects need not to be charged [HE]- Under 
rather general conditions (e. g., for non-magnetic objects 
in vacuum), the Casimir potential energy does not al- 
low for stable equilibrium positions of the interacting ob- 
jects [4]. This can be a practical disadvantage in systems 
where external (non-fluctuation) forces cannot be applied 
or fine-tuned to establish stability, especially in dynamic 
systems where the distance and hence the Casimir force 
changes in time. Nano-mechanical devices with closely 
spaced components fall into this class of systems [5] . Re- 
pulsive Casimir forces are known to exist if the space 
between the objects is filled by a dielectric with suitable 
contrast [3], but this is impractical in many situations. 

Repulsion can occur also in response to the prepara- 
tion of particles in distinct internal states, e.g. by optical 
excitation or coupling to heat baths of different temper- 
atures. In particular, Cohen and Mukamel |S] predicted 
a non-equilibrium repulsive force between molecules gen- 
erated by suitable detuning of the resonance frequencies. 
This suggests that for macroscopic objects held at dif- 
ferent temperatures similar repulsive effects should exist 
close to material resonances in the thermal region. How- 
ever, for macroscopic condensates the dielectric contrast 
is usually strong and non-equilibrium effects are compar- 
atively less significant than the equilibrium attraction. 
One should thus focus on sufficiently optically diluted 
materials to generate resonant Casimir repulsive forces 
between macroscopic bodies. The general formalism for 
dealing with non-equilibrium fluctuation effects between 
macroscopic bodies has been developed recently in [7HTT] . 

The aim of the present work is to explore theoreti- 
cally if and to what extent the above expectation can 
be realized and whether it can lead to stable equilibrium 



positions. We consider two dielectric slabs held at dif- 
ferent temperatures, and computed the pressure between 
them using the non-equilibrium extension of the Casimir- 
Lifshitz formula [7]. Employing a Lorentz-Drude dielec- 
tric response, we find that resonances and optical dilution 
indeed amplify repulsion sufficiently so that the total in- 
teraction can become repulsive, and equilibrium positions 
exist. We show that aerogels could be used to realize res- 
onant repulsion in practice since porosity can be used to 
tune reflectivity and resonances. The use of aerogels was 
indeed previously proposed to reduce the Casimir force 
[12], but not to generate repulsion. 

We consider two infinite parallel planar slabs S a , a = 
1,2, each consisting of a non-magnetic dielectric layer of 
thickness 6 and permittivity e a (uj), deposited on a thick 
glass substrate of dielectric permittivity e su b(cj). The 
slabs are held at temperatures T\ and T 2 , separated by 
a (vacuum) gap of width a. The Casimir pressure acting 
on the inside faces of the plates is given by [7] 
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where a is the Stcfan-Boltzmann constant. The last term 
in this equation is simply the classical thermal radia- 
tion result, which is independent of distance and mate- 
rial properties. In Eq. 0, P cq (T 1: T 2 ,a) = [P cq (T l7 a) + 
Pcq(T2,a)]/2 denotes the average of the equilibrium 
Casimir pressures at T\ and T 2 ; with P eq {T,a) given by 
the Lifshitz formula [14] . which for brevity is not repro- 
duced here. The novel non-equilibrium contributions are 
captured by the term AP noq (Xi, T2,a), which vanishes for 
T± = T 2 (it also contains a distance independent part). 
Upon decomposing AP ncq = AP™ + AP 1 ^ into prop- 
agating waves (PW) and evanescent waves (EW), one 



finds 
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where n(T) = [exp(b/i B T) - l]" 1 , fc z = ^w 2 /c 2 - k]_, 
Dp = 1 - r£V^ exp(2ifc z a), and P = M,N for the two 



polarizations. The reflection coefficients r P (w, k±) are 

(4) 



given by the well-known formulae for a two-layer slab 
„(«) __ r P (l, eg) + r P (e a , 6 sub ) e 2t ^° 



p l+r P (l,e Q )r P (e Q ,e sub )e M ?" ' 

where rp(e a ,ej,) are the Fresnel reflection coefficients 

e b (uj)q a (uj,k ± )-e a (uj)q b (uj,k ± ) 
e a (w)g b (u;, kj_) + e b (w)g a (u;, k±_ 



where q a = ^e a {u)uj 2 /c 2 - k\, and r M (e a ,eb) is ob- 
tained by replacing the e a and e b in Eq. ([5]) by one (but 
not in q a ). There is also an external pressure acting on 
the outside face of each plate which depends on the reflec- 
tivity of this face and the temperature of the environment 
T onv . In our numerical computations we assume that the 
external face of the glass substrate has been blackened, 
in which case the total pressure on plate a is given by 



P^(T U T 2 ,T C1 



la 



) = Pne q (Ti,T 2 , a) - — (T 4 + T 4 nv ). 

It is worth emphasizing that for T\ ^ T 2 , 
PW(T x ,T 2 ,T env ,a) ^ pWCTx.Ta.T^.o). Also, in equi- 
librium all distance independent terms vanish. 

The structure of the non-equilibrium force in Eq. M 
suggests that repulsion should exist out of thermal equi- 
librium. While the first term P eq (Ti,T 2 ,a) in Eq. Q 
is bound to be attractive, this is not so for the second 
term AP neq (Ti, T 2 ,a). As it can be seen from Eqs. (pi) 
and (pi , the quantity AP ncq (Ti , T 2 , a) changes sign if the 
temperatures T\ and T 2 are exchanged, and therefore its 
sign can be reversed by simply switching the tempera- 
tures of the plates. One notes also that AP neq (Ti, T 2l a) 



(i) 



<-> r 



(2) 



and 



is antisymmetric under the exchange r P 

vanishes for r p = r p . Therefore in order to take ad- 
vantage of this term to control the sign of the Casimir 
force it is mandatory to consider plates made of differ- 
ent materials. For real materials, both P cq (Ti, T 2 , a) and 
AP neq (Ti,T2,a) diverge as a~ 3 if a — > 0. In the follow- 
ing, we show that the sign of this asymptotic behavior 
can be made repulsive in certain cases. 



Motivated by the resonance-induced repulsion for mi- 
croscopic particles, we consider electric permittivities 
e a (oj) of the dielectric layers described by a Lorentz- 
Drude type model, as 



e a (w) 



C a w 2 a 



+ 



D a n 2 a 
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(An analogous two-oscillator model was used also for the 
permittivity e su b of the glass substrate, with the param- 
eters quoted in [2].) The first oscillator term (~ C a ) 
describes low-lying excitations of the materials; such low- 
lying polariton excitations in numerous dielectrics ac- 
count for sharp peaks in their dielectric functions in the 
far infrared region. Typical values for the resonance and 
relaxation frequencies are u> a = 10 13 — 10 14 rad/sec and 
7 Q = 10 11 — 10 12 rad/sec ;15j. The second oscillator 
term in Eq. ufh, proportional to D a , describes the con- 
tribution of core electrons. Excitation energies of core 
electrons are much larger, and characteristic values of 
il a are in the range 10 15 — 10 16 rad/sec. At and around 
room temperature core electrons are not thermally ex- 
cited (for T — 300K the characteristic thermal frequency 
ll>t = ksT/h is 3.9 xlO 13 rad/sec) and therefore their 
contribution to the thermal Casimir force AP neq is very 
small. However core electrons are important, as they 
strongly contribute to the average equilibrium Casimir 
force P eq , especially at sub-micron separations. 

Since the thermal Casimir force between two macro- 
scopic slabs should reduce in the dilute limit to the pair- 
wise interaction between their atoms, one expects that 
the resonant phenomena reported in Ref. [5] should be 
recovered if the material of the plates is sufficiently di- 
luted. In order to determine how large a dilution is nec- 
essary for this to happen, we investigated the behavior 
of the Casimir force under a rescaling of the amplitudes 
of the resonance peaks, i.e. C a — > C a /r, D a — > D a /r by 
a overall optical dilution parameter r > 1. 

We next report on numerical results based upon the 
above model for dielectrics, in which we set wi = 10 13 , 
7i = 72 = 10 11 , fti = ft 2 = 10 16 , and F x = T 2 = 5 x 10 14 
(all in rad/sec). In Fig. n] we plot the non-equilibrium 

(2) ~ 

normalized Casimir pressure Pnoq/P C q on slab 2 versus 
the ratio of the resonance frequencies uj 2 /tui, for a = 300 
nm, S = 5/im, C\ = 3, C 2 — 1.5, D\ — 1, and D 2 = 0.5. 
(Negative values of Pricq/Pcq correspond to repulsion.) 
For these parameters the non-equilibrium force is domi- 
nated by evanescent waves, whose skin depth is compara- 
ble to the separation a. This implies that for separations 
a <C S the force is practically independent of 6, but for 
a > 5 the features of the substrate influence the mag- 
nitude of the force significantly. The red dashed curves 
in Fig. [I] arc for T x = T cnv = 300K and T 2 = 600K, 
while blue solid curves are for T% = T cnv = 300K and 
T 2 = 150K. Three values of the dilution parameter t are 
displayed: r = 1 (x), r = 10 (+) and t = 20 (*). We 
see that for sufficiently high dilution, the Casimir force is 
strongly dependent on the ratio u 2 /w\ of the resonance 





FIG. 1: The normalized non-equilibrium Casimir pressure 
Pneq/Peq on slab 2 as a function of the ratio of the resonance 
frequencies W2/W1, for a = 300 nm, S = 5/xm, C\ = 3, C2 = 
1.5, D\ = 1, and D2 = 0.5. The red dashed curves are for 
Ti = T onv = 300K and T 2 = 600K, while blue solid curves are 
for Ti = T cnv = 300K and T 2 = 150K. Three values of the 
dilution parameter r are displayed: r = 1 (x), t = 10 (+) 
and r = 20 (*). Negative values of Pncq/Pcq correspond to 
repulsion. 



FIG. 2: The non-equilibrium normalized Casimir pressure 
Pneq/Peq on slab 2 as a function of plate separation (in mi- 



crons) for S - 5/um, Ti = T env = 300K, T 2 = 600K, r = 1 
and W2/W1 = 1.1 (x), r = 10 and W2/W1 = 1.05 (+), for 
r = 20 and W2/W1 = 1-03 (*). All other parameters are same 
as in Fig. [l] Dashed curves do not include the distance inde- 
pendent part of the pressure. Stable equilibrium positions are 
marked as SEP, unstable ones as UEP. The curve for r = 20 
approaches a negative values a — > 0. 



frequencies of the plates, displaying features analogous 
to those reported in Ref. [5] for the van der Waals inter- 
action between two 'atoms' coupled to baths at different 
temperatures. In particular, the Casimir force becomes 
repulsive if this ratio is suitably tuned to a value close to 
unity. We also find that the equilibrium force is rather 
insensitive to the tuning of resonances, and that the total 
force hardly depends on fi 1 /J7 2 - 

In Fig. [2] we plot the dependence of the non-equilibrium 
— . (2) ~ 

normalized Casimir pressure -Pn C q/-P e q on slab 2, as a 

function of plate separation a, at T\ = T cnv = 300K, 
T 2 = 600K, for r = 1 and w 2 M = 1.1 (x); r = 10 and 
w 2 /wi = 1.05 (+); t = 20 and w 2 /wi = 1.04 (*), all other 
parameters being same as in Fig. [I] The dashed curves 
do not include the distance independent part of the pres- 
sure, and are included to indicate that this component 
also changes sign upon dilution. We see that without 
dilution (r = 1) the non-equilibrium Casimir force is at- 
tractive for all displayed separations. By contrast, for 
large enough dilutions and for suitable values of ui^/uii, 
the force becomes repulsive in a wide range of separa- 
tions. The curve for t = 10 exhibits two points of zero 
force, one at a = 15 nm and another for a = 4.1 /xm, 
corresponding to an unstable equilibrium point (UEP) 
and a stable equilibrium point (SEP), respectively. For 
t = 20 there is a single turning point corresponding to 
a SEP at a = 3.3 /Ltm. There is no UEP point in this 

(2) ~ 

case as the ratio Pnoq/P B q approaches a negative value 
in the limit a — > 0, signifying that repulsion persists for 
arbitrarily small plate separations less than the SEP. 

As a means of achieving the dilution levels required to 
observe the resonance phenomena described above, we 
consider aerogels: highly porous materials fabricated by 
sol-gel techniques, starting from a variety of materials 



such as Si02, carbon, AI2O3, platinum etc. Aerogels with 
levels of porosity exceeding 99% can be realized nowa- 
days [IB]. In order to study the Casimir force between 
two aerogel plates we need an expression for the dielec- 
tric function e(w) of the aerogel, valid in the wide range 
of frequencies relevant for the Casimir effect. For sepa- 
rations a larger than the pore size (typically of the order 
of 100 nm or less), an effective medium approach can be 
used in which the aerogel permittivity e(u>) is obtained 
from the Maxwell- Garnett equation [IT] as 



1 



(1 



1 



(8) 



Here, e(w) is the permittivity of the solid fraction of the 
aerogel, and < < 1 is the porosity. Equation ffity is 
justified if the solid fraction is well separated by the host 
material (air), i.e., when (f> is sufficiently close to one. 
Using again a model of the form in Eq. (TtJ) for e a (ui) , one 
finds that e a (ui) has a resonance at the frequency 



1 



(p a C a 
t> a D a + 



1/2 



(9) 



where we assumed w Q -C Q a . According to Eq. J9J, the 
frequency of the aerogel resonance is blue-shifted with 
respect to u> a , and as <j) a is varied from zero to one, <2> a 
sweeps the range from w a to [1 + C a /(D a + S)] 1 / 2 ^. 
The dependence of the resonance frequency w Q on the 
porosity is welcome, because it gives us the possibility 
of tuning the ratio of the resonance frequencies for the 
two plates by simply choosing appropriate values for the 
porosities <fii and (\>i of the plates. 

As an example, we consider two aerogel layers of thick- 
ness 6 = 5/xm, deposited on a thick glass substrate, with a 
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FIG. 3: The pressure ratio Pn e q/-P Gq for two aerogel layers 
(<5 = 5 /im) on a glass substrate as a function of the porosity 
4>\ of plate one. The porosity fa of plate two is 0.95 (+), 0.9 
(x) and 0.8 (*). The dielectric functions e a (uu) of the host 
materials are as in Eq. |m|, with C\ — 1, C2 = 3, D% — D2 = 
0.5, W2/k>i=0.84. All curves are for a separation of 200 nm 
and for Ti = T cnv = 300 K, T 2 = 600 K. 




FIG. 4: The pressure ratio Pneq/Pcq for two aerogel layers 
(<5 = 5 /im) on a glass substrate as a function of plate separa- 
tion (in microns), for fa = 0.95 and ^2=0. 95 (+), 0.9 (x) and 
0.8 (*). All parameters for the aerogel plates are same as in 
Fig. [3] The dashed lines do not include the separation inde- 
pendent part of the pressure. The stable equilibrium position 
is marked by SEP. 



blackened outer surface. The dielectric functions e a (u>) of 
the host materials are as in Eq. Q, with C\ = 1, C2 = 3, 
D x = D 2 = 0.5, w 2 /Vi=0.84. In Fig. pj] we plot the ra- 

tio Pncq/Peq as a function of the porosity fa of the first 
plate. The porosity fa of the second plate is 0.95 (+), 0.9 
(x) and 0.8 (*). All curves in Fig. [3] are for a separation 
of 200 nm, at T x = T cnv = 300 K and T 2 = 600 K. The 
force can indeed be made repulsive by suitably adjusting 
the porosities of the two plates. In absolute terms, the 
Casimir force is small; e.g. for fa — 0.77 and fa = 0.8 



we find P, 



(2) _ 



0.55 x 10" J Pa. 



In Fig. UJwe plot dependence of the ratio Pn Cq /P C q on 
the separation (in microns), for fa = 0.95 and </> 2 =0.95 
(+), 02=0.9 (x) and 2 =O.8 (*). All parameters for the 
aerogel plates are the same as in Fig. [3] The dashed 
lines do not include the separation independent part of 
the pressure. We find that for fa — fa — 0.95 a stable 
equilibrium point exists (marked as SEP in Fig. Rj) . 

While fluctuation induced forces are generally attrac- 
tive, repulsive forces can be obtained between atoms pre- 
pared in different excited states. Coupling of atoms to 
thermal baths at different temperatures can in princi- 
ple produce population of states that lead to repulsion. 
We show that a similar non-equilibrium repulsive contri- 
bution to pressure also arises from the interplay of res- 
onances for two macroscopic slabs held at appropriate 
distinct temperatures. However, to observe a net repul- 
sion between slabs one must overcome an ever present 
attractive force arising from the dielectric contrast of the 
condensed bodies from the intervening vacuum. The lat- 
ter can be reduced by dilution, and we have shown that 
aerogels provide a material where a net repulsion at sub- 
micron separations can be achieved, leading to a stable 
equilibrium point. 
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